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Deterministic generation of arbitrary photonic states assisted by dissipation. 
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A scheme to utilize atom-like emitters coupled to nanophotonic waveguides is proposed for the generation 
of many-body entangled states and for the reversible mapping of these states of matter to photonic states of an 
optical pulse in the waveguide. Our protocol makes use of decoherence-free subspaces (DFS) for the atomic 
emitters with coherent evolution within the DFS enforced by strong dissipative coupling to the waveguide. By 
switching from subradiant to superradiant states, entangled atomic states are mapped to photonic states with 
high fidelity. An implementation using ultracold atoms coupled to a photonic crystal waveguide is discussed. 

PACS numbers: 42.50.-p,03.67.Bg,42.50.Ex 


Introduction. Recent experimental and theoretical work on 
optical emitters coupled to one-dimensional waveguides has 
opened new avenues to investigate light-matter interactions 
[1-19]. Particularly promising are the setups where atoms 
are strongly coupled to structured dielectrics [6-10], where 
large Purcell factors have been predicted [20, 21]. Further¬ 
more, collective effects can be enhanced by placing the atoms 
at particular positions [15, 16, 22-28]. The combination of 
atom-like emitters and nanophotonic waveguides may enable 
new regimes for the interaction of light and matter, leading 
to technologies that outperform current ones and qualitatively 
different physics. In this work we investigate the possibility of 
using atom nanophotonics interfaces to tailor arbitrary states 
for propagating photons on demand. We predict large fideli¬ 
ties even for relatively large numbers of photons, something 
which has been impossible to achieve with other platforms in 
the optical domain. 

Our proposal uses three-level systems (with levels 
{|g), |v), 1^)}), where one of the optical transitions (|g) \e)) 

is strongly coupled to a one dimensional (ID) waveguide (see 
Fig.l(a-b)). We denote by Pid the Purcell factor correspond¬ 
ing to that transition; i.e., the ratio of the emission rate into 
the particular waveguide mode. Fid, and the one for all other 
modes, F*. The atoms will be separated by distances propor¬ 
tional to the optical wavelength = 2n jqifof), where q{(B) 
is determined by the waveguide dispersion relation. Depend¬ 
ing on their internal state, atoms may experience a collective 
decay into the waveguide, or become completely decoupled 
from it. The latter occurs if they are in a decoherence free 
subspace (DFS) [29, 30], that is, if they are the states from the 
nullspace of the collective decay of the transition \e) -> \g). 
Our protocol consists of two steps: in the first one, we gener¬ 
ate certain decoherence-free states, j'P/)), by driving the atoms 
with lasers and using the collective quantum Zeno effect [31- 
33] within the DFS with a fidelity 1 —/q oc m/y/Pyo, where 
m is the maximum number of photons we want to generate; 
in the second one, a laser pulse takes the atomic state out 
of the DFS so that atoms collectively emit into the waveg¬ 
uide, creating the desired state of a single propagating mode. 
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Figure 1. (a) Scheme of the proposal: N atoms equally spaced, 
Zn = wAa = nlizlq^,, with n eN, plus 1 ancilla atom coupled to a ID 
photonic bath that can be separated by a larger distance as long as it 
satisfies Zn = nX^. The ancilla atom must be individually addressed, 
(b) Atomic A-scheme considered: the dipole transition \g) ^ \e) 
is coupled to the reservoir modes, aq. A Raman laser controls the 
transition \s) ^ \e) with intensity, kin and detuning A„. A second 
field controls \s) ^ |g), with frequency cOg. (c) Relevant states of the 
atomic ensemble for our three step protocol: (I) generation of super¬ 
positions of symmetric states, |D^), in the two metastable states \g) 

and \s) by using the excited dark states (II) Flipping |^) ^ \e) 

to generate the superradiant state \Sm), which decays rapidly (III) to 
the desired entangled photonic state. 


As we show, the second process has an (in)fidelity 
I—f 2 oc /NPyo, where N is the number of atoms. 

Coupling to a common environment The atom-photon 
hamiltonian of these systems is given by H = Hq ^ Hi, 
where Hq = F/qb + Afield is the free energy term, with F/qb = 
Lnll ^nd //field = Y^gCOquloq, (using H = 1) 

where cOa is the two-level system energy, (Oq is the field dis¬ 
persion relation of the ID photonic modes, and = \i)n{j\- 
We consider only the coupling to a single polarization as jus¬ 
tified for suitable dielectric waveguide modes [20], that is, a 
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dipolar coupling of the form 

= ( 1 ) 

n,q 

with gq the single-photon coupling constant to the mode 
of interest and where we have used the rotating wave ap¬ 
proximation. Under the condition that the ID-bath degrees 
of freedom have a much faster relaxation timescale than 
for the atomic system, the atoms are described by a den¬ 
sity matrix, p, which in the Born-Markov limit, is governed 
by a master equation [22, 23, 34] of the form: dpjdt = 
(a«,pa,^-pa,^a«,) +H.c., where Fid 
is the renormalized spontaneous emission rate due to the in¬ 
teraction with the ID photonic reservoir. 

Emergence of subradiant and superradiant states. By ap¬ 
propriately choosing the atomic positions, e.g., Zn = = 

nlnjqicof), with ^ G N, the coherent atom-atom interactions 
are eliminated [35] and the effective interaction yields a pure 
Dicke model [36] decay described by 

^d(P) = ^ - SegSgeP) + H.C., (2) 

where we have introduced the collective spin operators Sij = 
ofj. One of the assets of the model described by Eq. 2 
is the emergence of sub and supervdidmni states. The excited 
states with m atoms in |^) that are symmetric under the per¬ 
mutation of atoms, denoted by \Sm), are superradiant with an 
enhanced decay rate proportional (at least) to the atom num¬ 
ber N, and are unique for each m. On the other hand, the 
states satisfying = 0 are dark states of the Liouvillian 

of Eq. 2, and are therefore decoupled from collective dissipa¬ 
tion. These dark states span the DES that is highly degenerate 
for m > 1. 

The atomic entangled states that have to be created in the 
first step of our protocol are very peculiar, as: i) they have to 
belong to the DES; ii) they must be prepared within that sub¬ 
space to avoid dissipation; iii) and they must have a special 
form such that they can be easily mapped in the second step 
to the appropriate superradiant states, in order to generate ar¬ 
bitrary superpositions of the photonic states in the waveguide. 
Our strategy consists of first, identifying states, denoted by 
\Dm), in the subspace spanned by the ground levels g, and 5-, 
which can be mapped one-to-one to a basis \Sm) of superra¬ 
diant states using a simple laser pulse, and which in turn give 
rise to m photons in the waveguide via superradiance. Then, 
we use a more sophisticated scheme within the DES to gener¬ 
ate superpositions of |D^), which requires only m steps. 

By introducing an auxiliary metastable state, \s), as de¬ 
picted in Eig.l(b), the candidate states to map to superradi¬ 
ant ones are the symmetric Dicke states |D^) oc sym{|5')'^'^ (g) 
as they can be turned superradiant by switching 
\s) \e). Having identified the target, |D^), we need to find 

ways to build efficiently arbitrary superpositions. One gen¬ 
eral strategy already explored in the literature consists of im¬ 
plementing one and two-qubit universal gates within the DES 


[29-32], which in principle allows to build any arbitrary su¬ 
perposition. However, the number of steps increases rapidly 
with the number of atoms. Another possibility is to use su¬ 
perposition of Zeeman states of 1 atom and rapid adiabatic 
passage methods [37-39], but it is limited by the number of 
metastable states available. Eor small system sizes other adi¬ 
abatic passage methods have been numerically studied [40]; 
however, they do not work for superposition states. In our 
proposal, we use the collective character of the interaction to 
deterministically generate arbitrary A-qubit states for which 
the number of steps is independent of the number of atoms. 

Controlling atomic states under strong dissipation. 

The scheme that we use is depicted in Eig. 1(a): we consider 
a system of A + 1 emitters, in which we aim to generate \Dm) 
in the first N emitters using the ancilla as an auxiliary state. 
As the \Dm) ’s are invariant under the permutation of the first N 
atoms, we choose the control fields with the same symmetry. 
They can be described by the following hamiltonian (in the 
interaction picture with respect to iTqb): 



(3) 


^las — 


\ A^+1 

H.c. +A, £ (7" , 

/ n=\ 


where allows to control both the emitter state and the 
coupling to the ID-reservoir, while allows to control 
the atomic states of the ancilla independently of the cou¬ 
pling to the ID-reservoir. We are interested in working in 
the regime of strong collective dissipation, where AEid ^ 
^r,^anc 5 ^c,Ag- In this situation, the ID-bath is continu¬ 
ously monitoring the collective atomic state, as in the Quan¬ 
tum Zeno regime [29-32], and projecting the atomic state into 
the DES of the Liouvillian ^d- Eormally, we obtain the effec¬ 
tive dynamics within the DES by treating the control fields 
and dissipation into other decay channels as a perturbation to 
the collective dissipation [41]. Eor the ideal protocol, we 
consider the first order, that is the action of the control fields 
projected into the DES. The errors due to higher order terms 
are discussed later. 

Preparation ofmany-body entangled states: ideal protocol. 
Due to the symmetry of the problem, it is convenient to intro¬ 
duce the following notation to describe any symmetric state 
over A atoms: 


\F,n,k) - sym{|.)®'" ® ® , (4) 


that embeds both |D^) = \Fmp) and \Sm) = |Tb,m)- In general, 
the DES of the Liouvillian grows exponentially with the 
number of atoms, but for each m, only three of these states 
fulfill the permutation symmetry of our system [41]. These 
states are: 


=|^m-1.0) 0 k)A , I'I'f = \Fn,fi) 0 |g)A , 


(5) 




0 k>A - 0 Is)*. 
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Figure 2. Scheme for the preparation of arbitrary superpositions of 
symmetric Dicke states |D^) of N atoms. Alternating between Gx 

gates on the ancilla and two-photon transitions via builds up 

excitations step by step (from left to right). To make the ladder-like 
structure clearer, we added a term for the ancilla in the energy scale. 


where \\I/)a denotes the state of the ancilla atom and Nm = 
N — In Fig. 2 we sketch the protocol steps by a di¬ 

agram of the projected hamiltonian into the DFS. It consists 
of two parts, which are applied mmax times to reach any su¬ 
perposition of states over the first N atoms j'F/)) (g) Ig)^ = 
(Lm=0i ^m\Dm)) ^ |^)a with maximum mmax excitations from 
the initial state 

1. Use He to flip the ancilla state |g)A ^ |‘S')a- This transi¬ 
tion makes: ^ 

2. Two-photon transition ^ It can be 

shown that the dark states corresponding to a given ex¬ 
citation m form an effective A-scheme within the DFS 
via the far detuned state [41]. The two-photon 

transition can be made on-resonance if the Rabi cou¬ 
plings are chosen such that |^2anc| = which 

is possible because we demanded individual addressing 
of the ancilla. Notice that if this condition is not im¬ 
posed, the experience different Stark-shifts that 

ultimately spoil the two-photon process. The effective 
Hamiltonian for the m-th excitation is 

//D = QW/2|»p('”)){'pf)|+H.c., (6) 

where and where the Stark shifts 

compensate each other and can thus be neglected. 

As depicted in Fig. 2, the combination of both steps gives 
rise to a ladder-like structure, which can be used to build any 
arbitrary superpositions state of |D^) from the ground state 
One way of finding the pulse sequence to create any 
arbitrary superposition consists of removing the excitations 
in the system, starting from step by step until only 

is populated [33]. Remarkably, the number of steps 
required depends only linearly on the number of excitations 
^max, independently of the number of atoms, in stark contrast 
to building the states out of one and two-qubit gates [31]. 


Shaping an arbitrary field of light using strong dissipation. 
Apart from the intrinsic interest of generating arbitrary super¬ 
positions I'R/)), another advantage is that our protocol can be 
used to generate a superposition of bosonic states in the ID- 
bath by dissipative means. To make sure the ancilla atom can 
be neglected, we flip the ancilla state |g)A ^ |^)a and apply 
no fields to it. In order to map to the superradiant state of N 
atoms, we apply a fast resonant ;r-pulse (A^ = 0 and ^ 
A^Fid) on the N emitters to switch all \s)n ^ \e)n, thus gen¬ 
erating the superposition of \Sm) ^ sym{|^)®"^ (g) \g)^^~^}. 
Due to their superradiant character, the \Sm) decay completely 
to ID-reservoir modes. Because Hi conserves the number of 
excitations, the superradiant state of m-excitations, decays to 
Fock-state of m-photons [41]: 

\Sm) ^ Im^) = £ ’ (7) 

{q} 


where {q} = {q \,..., is the set of relevant momenta of the 
m-photon Fock state, where each of them run over the whole 
Brillouin Zone qt G B.Z.. The scattering amplitude 
can be calculated using a generalized input-output formalism 
[41-43] and quantum regression theorem [44]: 


m 


=n 

r=l 


\g^JrNr e 

i(52/=l ^qi ~ + rY\Y)Nr/7. 


+P[M] (8) 


for sufficiently large times t ^ I/A^^Fid when the atomic 
state has decayed completely and defining P[{^}] as all the 
permutations of {q}. Notice that the only dependence on t in 
this case enters through: e' -ilr=l , which describes the cen¬ 
ter of mass motion of the wavepacket when going to the real 
space. In the low excitation regime, one can either use the 
Holstein-Primakoff approximation [45] or directly substitute 
An -G A in the expression above. In both cases, we arrive to 
[41]: 

yi(%-ffla)+riD^/2’ 



that has a Lorentzian shape centered at (Oa with bandwidth 
FidA/ 2. Substituting -g into the definition of 

yields a linear Fock state denoted by In prin¬ 

ciple, the emission into the waveguide is bidirectional (±^), 
but combining both fields in phase, e.g., by a placing a mirror 
at an appropiate distance, or by engineering the atom-photon 
coupling appropriately [3, 4, 14], it is possible to achieve 
emission in one-direction only. Furthermore, by shaping the 
pulse, Ylx{t) (within a bandwidth < AFid) we generate any 
desired shape of the output photonic state, e.g., to create a 
time-symmetric photonic state [46] that ensures the reversibil¬ 
ity of the process when mapping the photonic state to another 
sample. Moreover, because of the linearity of the calculation 
of Aj^j with respect to the input state [41], superpositions of 
atomic states decay to superpositions: 


Wmax mmax 

I'Vd) ^ £ d^\s^) ^ I'Pb) = £ , (10) 

m=0 m=0 
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that will be generated in a single-mode wavepacket as long as 
N ^ m because: (m|^j = I—j (20A^^) + 0{m^ /N'^) 

[41]. 

Preparation ofmany-body entangled states: error analysis. 
So far, we have only discussed the ideal protocol without con¬ 
sidering detrimental effects, e.g., spontaneous emission into 
all other modes with rate F*. For the error in the prepara¬ 
tion of the many-body entangled state, we derive an error rate 
£ from perturbation theory, which, together with the time of 
the operation T, gives an approximation of the infidelity of 
the state, 1—Fi ^ T£, where the fidelity Fi is calculated with 
respect to the target state I'Fd), i.e., Fi = ('Fd|p(t)|'Fd). 

The dominant errors assuming N m and Fid ^ ^ F* 

(see Ref. [41] for complete expressions) come from: i) the 
spontaneously emitted photons from to decay channels 

other than the waveguide. As we are using an off-resonant 

Raman transition, this error scales as £\ ^ F* • h) Other 
errors appear due to the photons emitted from the small pop¬ 
ulations of superradiant states. We estimate the rate of these 
errors by taking into account the second order corrections of 
the projected hamiltonian [41], which are finally given by 

m\£l 1 ^ 

£2 ~ NT ID 4(a|+(atid)^) ’ ^P’ 1^^ infidelity for the 

m-th step of the process, which takes T = ~ 

full population transfer, is 


l-F, 


im) 


2 



( 11 ) 


that is optimised for A^^opt = \/r*FiD, which yields a scaling: 

1 — opt ^ 1 /\/Fid- Interestingly, the scaling depends neither 
on the number of atoms, N, nor the number of excitations, m. 
In order to create any superposition I'Fd), we require mmax 
steps, thus, the total infidelity of the first part of the protocol 
is 


. ^ ^max 

vFid 

that can be improved via post-selection of the states condition 
on detecting no photons to the waveguide. [47] 

Preparation of many-body entangled states: numerical 
analysis. To validate the scaling analysis, we study numeri¬ 
cally (without any approximation) the preparation of two rel¬ 
evant sets of states: i) the general class of states |F>m); ii) 
the superpositions ^ (|Z)o) + |F>^)). First, recall 

that the symmetry conditions that we found from our anal¬ 
ysis of the ideal situation tell us that the relevant Hilbert 
space does not depend on the atom number, N, but only 
on the maximum number of excitations, mmax [41]. The 
atom number only enters on the two-photon resonance con¬ 
dition, that fixes Using this restricted Hilbert space, 

we use a non-hermitian evolution determined by the coher¬ 
ent driving plus the imaginary energies determined by the 
coupling to the waveguide through the Sge operator, to give 
F/eff = Ffias — iFiDAegAge/2 — iF*5'ee/2. 



100 200 500 1000 2000 5000 
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Figure 3. (a) [(b)] Infidelities as a function of Pid for generating 
\Dm) [|<l>m) = Tfd^o) + \Dm))] up to mmax = 5 excitations. The 
dots correspond to the numerical fidelities, whereas the solid lines 
depict the 1 / v^Pid scaling from Eq. 12. 


For the generation of Fock states, the pulse sequence con¬ 
sist of a complete transfer of populations in each step of Fig. 2: 

^ ^-^ which can be ensured by fix¬ 
ing the time of interaction, t, to tTlc = Tt = n) for the 

microwave (Raman) transitions. For a more complicated state, 
such as the |<I>^) the pulse sequence can be calculated numeri¬ 
cally. In Fig. 3(a-b), we show the numerical fidelities obtained 
when fixing the off-resonant transition to the optimal A^ opt 
that we found from our error analysis. The plots show how 
our general arguments give the correct scaling oc 1 /^P^d- 
Generation of Photonic States: Error Analysis. The errors 
from deviations from the Holstein-Primakoff approximation, 
which leads to Eq. 9, have already been discussed. Here, we 
take into account the emission into other modes through F* 
to obtain an estimate for the fidelity, that for a state with m- 
excitations is given by mF*. Then, the infidelity of an op¬ 
eration with average time is approximately 1 — F^^"^ = 

. When there are mmax excitations in the system, one ob¬ 
tains an infidelity for the mapping to the photonic state of the 
process: 


^max^* _ ^max ^max 

AFid 


(13) 


Remarkably, the dissipative character of this mapping al¬ 
lows for the deterministic and efficient generation of (arbi¬ 
trary superpositions of) photonic states, e.g., Fock states, that 
typically must be generated probabilistically [48, 49] with an 
exponential scaling with the number of excitations or via non¬ 
linear interactions [33]. 

Implementation. In dielectric photonic crystal waveguides, 
Fid/F a = ^ngOj{lAm), where ng = c/vg the group in- 
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dex, a = 1(271) the radiative cross-section, the ef¬ 

fective mode area, Fa the vaccuum emission rate and ^ an 
adimensional factor coming from cavity enhancement. Cur¬ 
rent values for Cs atoms (Aq = 894 nm, ra/2;r = 5.02 MHz) 
and SiN alligator waveguides [6, 10] have A^ ^ 0.2 /rm^, 

^ 10, ^ 5 and Q-factors of 10^. Depending on the re¬ 

duction of spontaneous emission, F* = aFa, these numbers 
lead to PiD ~ 50/a. Intrinsic losses in the dielectric (finite 
2-factors) and reduced set finite propagation lengths of 
waveguide modes, Fprop/Aa ~ Ql{2Kng), which is > 10^ for 
state-of-art SiN values [6, 10]. Retardation effects also set 
a maximum number of atoms for which superradiant atom- 
photon mapping could be observed, e.g., assuming a separa¬ 
tion Aa/2, then A^Fid < 2vg/(A^Aa), which for current struc¬ 
tures leads to < 500 atoms. Possible ways of overcom¬ 
ing retardation are to increase Fid not only by Vg but through 
cavity enhancement, e.g., by using mirrors [10]; or directly, 
by doing the atom-photon mapping off-resonantly, which de¬ 
creases effectively Fid (and F* by the same amount) relax¬ 
ing the retardation requirements. Other potential problems 
that we have neglected include i) imperfect atomic separa¬ 
tions limited ultimately by quantum center of mass wavepack- 
ets and atomic motion and ii) group velocity dispersion that 
distorts the wavepacket as they propagate through the waveg¬ 
uide. These effects should be studied in detail to evaluate and 
minimize the impact on the performance on the protocols de¬ 
scribed. 

Conclusions. We have presented a protocol to generate 
deterministic superpositions of many-body entangled atomic 
states in the presence of strong dissipation. Remarkably, the 
infidelities in the preparation of complex superposition states 
increase only linearly with the number of excitations of the 
system (m), inversely with the square root of the Purcell Fac¬ 
tor (1 /v^Pid), and are independent of the number of atoms 
considered. Furthermore, we have shown how to map these 
atomic states to photonic states with a very efficient scaling 
that depends linearly on the inverse collective Purcell Factor 
1 / (A^Pid) and how to engineer a time-symmetric wavepacket 
that guarantees the reversibility of the mapping. A possi¬ 
ble outlook of this work is to use other metastable states and 
different polarizations to generate entangled photonic NOON 
states in directions [3, 4, 14] which is very relevant for 
quantum metrology applications [49, 50]. 

Acknowledgements - We gratefully acknowledge discus¬ 
sions with T. Shi. The work of AGT, VP and JIC was 
funded by the European Union integrated project Simulators 
and Interfaces with Quantum Systems (SIQS). AGT also ac¬ 
knowledges support from Alexander Von Humboldt Foun¬ 
dation and Intra-European Eellowship NanoQuIS (625955). 
JIC acknowledges support as a Moore Distinguished Scholar. 
DEC acknowledges support from Eundacio Privada Cellex 
Barcelona, Marie Curie CIG ATOMNANO, and ERC Starting 
Grant EoQAL. HJK acknowledges funding by the Institute of 
Quantum Information and Matter, a National Science Eunda- 
tion (NSE) Physics Erontier Center with support of the Moore 
Eoundation, by the Air Eorce Office of Scientific Research, 


Quantum Memories in Photon-Atomic-Solid State Systems 
(QuMPASS) Multidisciplinary University Research Initiative 
(MURI), by the Department of Defense National Security Sci¬ 
ence and Engineering Eaculty Eellows (DoD NSSEEE) pro¬ 
gram, by NSE PHY 1205729 and support as a Max Planck In¬ 
stitute for Quantum Optics Distinguished Scholar. 


6 


Supplemental Material: Deterministic generation of many-body entanglement and photonic states 

assisted by dissipation. 


PREPARATION OF |£)„). 
Coupling to a common environment. 


Let us consider first the hamiltonian part describing the atom-photon coupling for an optical transition \g) <—> \e), coupled 
to ID photonic reservoir. In this case, the hamiltonian is given hy H = Ho-\-Hi, where Hq is the free term: Hq = i^qb + Afield, 
(using h= 1) 

A^+l 

^qb = tOa ^ CJgg, Afield = (SMI) 

n=\ q 

where tOa is the atomic transition energy, CDq is the field frequency from the dispersion relation of the waveguide modes. We 
consider that we are coupled to a single polarization to focus on the most relevant physics of our work which can be justified for 
appropriately designed dielectric waveguides [20] . We consider a dipolar coupling of the form 

iV+l 

= E {^geE{Zn) + H.c.) , (SM2) 

n=l 

with E{z) = and gq the single-photon coupling constant. In the conditions where the ID-bath degrees of freedom 

have a much faster relaxation timescale than the system, we can describe our atomic system via a density matrix, p, which in 
Born-Markov limit, is governed by a master equation: dp/dt = ^d(p) [22, 23, 34], with the superoperator 

^d(P) = E i(y"eP<g -p<yX) + H.c., (SMS) 

n,m 


where Fn^rn = with Fid the renormalized space spontaneous decay rate of the atoms due to the interaction with 

the ID photonic reservoir. For completeness, it is interesting to write the connection between Fid and gq that can be easily 
computed from [34] : 


rio = 27tY^\gg\^5{0)^-0)q)=L [ dq\gq\^8{(i)^ - (O, 

^ J —oo 


Vg{(0^) 


(SM4) 


where we have introduced L the quantization length of the guided modes and used the fact that COq = CO-q and \gq\ = \g-q\- 


Emergence of subradiant and superradiant states. 

In the main manuscript we stipulated a homogeneous coupling to the environment, that can be achieved naturally with a ID 
reservoir by choosing appropriately the atomic positions, i.e., Zn = = n27tlq{ci)^), with n eN, along the waveguide. With 

this choice, the effective interaction induced by reservoir modes yields a pure Dicke model [36] decay described by 

^d(P) = ^ {SgepSeg - SegSgeP) + H.C., (SMS) 

where we have introduced the following notation for the collective spin operators Sij = Seg is just the collective 

operator for the spin dipole Geg. One of the assets of the model described by Eq. SMS is the emergence of suh and supervdidmni 
states as depicted in Fig. 1(c), which can be seen easily by examining Sge\J^mj^ aj) = Y^(/ + m/)(/ —m/ + l)|/,m/ — 1, a/) 
in the collective angular momentum basis aj)} with J = N/2^N/2 — 1, ... and mj = —and aj is an 

index that takes into account the degeneracy of the states (that we drop from here on as it does not play an important role 
for what we will describe). The states satisfying = 0 (i.e., mj = —J) are dark states of the Liovillian of Eq. SMS 

and they form a so-called Decoherence-Free Subspace (DFS)[29, 30], and are therefore uncoupled from dissipation, whereas 
the states with J = N/2^ mj > —N/2 are super-radiant with an enhanced decay rate proportional to the atom number N. It 
is interesting to emphasize that we can find similar physics, i.e., pure Dicke model, using an atomic configuration where the 
Zn = ^Aa/2 = n7tlq{(0^) in G N). The difference between the two configurations is the symmetry of the super/subradiant states, 
which has to be taken into account in the generation of these states. 
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Controlling atomic states under strong dissipation. 


Our first goal is to generate particular superpositions of atomic states within the DFS; therefore we need to include in our 
system some additional fields that allow to control the individual atomic states. We find convenient to introduce an extra 
auxiliary level \s)n, as sketched in Fig. 1(b) of the main text and use the following fields to control the atomic state: 


f Cl ^ Cl \ ^+1 

Hus ={ -yL<e+ ^ +h.c. + A, £ ct; 

\ ^ n=l ^ / n=l 

= Y<+'+h.c., 


(SM6) 

(SM7) 


where Ffias allows to control both the coupling between atom and ID-reservoir, while allows to control the atomic state of 
the ancilla atom independently of the coupling to the reservoir. We are interested in working in the regime of strong collective 
dissipation, where A^Fid ^ In this situation, we can intuitively consider that the ID bath is continuously 

monitoring the atomic state, as in the Quantum Zeno regime [29-32], and projecting the atomic state into the DFS of the 
Liouvillian Notice that when including the extra auxiliary level \s), the DFS = {j'F) : Sgej'F) = 0} contains all superposition 
states of the metastable states \g) and \s) plus all the excited states from the nullspace of the collective dipole Sge. Formally, we 
obtain effective dynamics within the DFS by using a projector operator P satisfying: P^d = = 0, and its orthogonal part: 

Q = 1 — P. Using these projectors, one can formally integrate out the fast dynamics outside of the DFS, described by Qp, and 
obtain effective dynamics of the atomic system within the DFS given by [34] 


1 t“3 

Pp = PWPp - PWQ^QWPp + 0[^]. (SMS) 

Afb rfo 

where W is any perturbation acting on the atomic system (with relevant time scale t), e.g., W = Ffias +F/c and tFid must be <C 1 
such that it is a good approximation. 


Preparation of many-body entangled states: neglecting losses. 

As explained in the main manuscript, the first step of our protocol consists of creating a certain class of states that must satisfy: 
i) they must be easily mapped to the superradiant states of the atomic ensemble; ii) they must be created using only states within 
the DFS such that we can avoid the dissipation induced by the waveguide. We propose in Fig. 1(a) of the main manuscript a 
configuration of N atoms and a separately addressable ancilla atom. Because of the high symmetry of the states that we aim to 
create among the first N atoms, it is convenient to introduce the following notation to describe any symmetric combination of 
states over these atoms: 


(SM9) 

where ^{m,k) = („,(. multinomial coefficient that gives the normalization of these states. This notation embeds 

the many-body entangled states that we aim to create, i.e., \Ffn^o) = |D^) and iFb^m) = \Sm)- For the ancilla atom, we use a 
notation \ \1/)a- 

In Eqs. SM6-SM7 we introduced the laser configuration that we use, namely, a symmetric excitation over the first N atoms 
and a different one for the ancilla. Interestingly, the combination of this configuration and the collective dissipation imposes 
certain symmetry conditions on the states that can exist within the DFS, namely, the atomic states with excited states (\e)) must 
be symmetric under any permutation of the first N atoms and antisymmetric with the permutation of the ancilla. This reduces the 
exponential Hilbert space of all states (3^) to a set of 5m relevant states, which are depicted in Fig. SMI. For example, among 
all the combination of states with one atomic excited state, for each m only one of them (denoted by belongs to the 

DFS. Moreover, the states of combinations of the ancilla with |F^,2) are all superradiant as none of them can fulfil the symmetry 
to be within the DFS of the collective dipole Seg. 

Now, let us explain how to generate the |D^) using the tools that we have introduced. As shown in Fig. SMl(b), for each 
m, there exists an effective A-scheme within the DFS of the Liouvillian that couples the states = |D^-i) ( 8 ) \s)a to 

= \Dm) 0 |g)A, through an excited state that also belongs to the DFS. If no projection into the DFS is considered. 
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(a) \Fm-Z,2) ^ \s)jFm-2,l) ^ |e)^|F^-2,2) 0 |^/^ 







cf/l 

1 


> 
d ^ 


d 

i) \9\ 

-2, 


■i,o) 0 |e)J^m-i 


> 




0^7 




Oc 


(b) 


|i^m-l,o) < 8 > 1^?) \Fm-l,o) 0 |s), \Fm,o) 0 \g)j^ \Fm,o) 0 |s)^ 


(X5e5^^^) OC Seg[^^) OC Seg\)^) 


ixi”>) I'l-r') ixr') 


Mv 






|g,(m-l)) 


I^M) 




|^("^+i)) 


Figure SMI. (a) Product of the Hilbert space of N permutation invariant emitters and of one ancilla for m excitations, that is m emitters in 
state \s) or \e). (b) Separation of whole Hilbert space (for m excitations) into DFS states (blue background) and non-DFS states, which makes 
obvious the emergence of the effective A-type transitions within the DFS. 


the states are coupled with excited states as follows: 




f^m-l 


,o) k)A, 


(SMIO) 


Interestingly, it is possible to write Eqs. SMIO separating the contributions of the states in and out of the DFS: 






(SMll) 




^ I^p(m) 


s/ynNfYi I (m)\ 

■M)\ 


where we have introduced the notation — m + 1, and is a state within the DFS that couples both and is 

given by: 

\X^^g) are two states outside the DFS defined as 

|xi”’^) = |F„_2,i)0|5)A, (SM12) 

and can be shown to have an enhanced decay rate of F^ = {Nm + l)riD, by looking at the action of the collective operator Seg 
on these states. We discuss first the effect considering only perturbations up to first order within the Zeno dynamics. Therefore, 
the super-radiant states can be neglected as they are only virtually populated due to their enhanced decay rate. Thus, we first 
consider the effective A system, with effective Raman intensities given by (see Fig. SMl(b)): 


^(m) ^ = _ 


(m) \ _ ^anc 


Nyn 


A^m + 1’ 


(SM13) 




2 V + l 

where we see the importance of addressing the ancilla atom individually from the other N emitters and keeping Q.r 7^ ^anc, 
as we can now set them such that by choosing: Inland = This choice allows us to compensate 

the different Stark-shifts that are introduced with the projection P and yields an off-resonant two-photon transition with Rabi 
frequency: 


|qW| = 


IQ, 


2Ag N,n +1 


(SM14) 
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Notice then, that by flipping the state of the ancilla with one also flips ^ which re-initializes the 

process. Thus, by using a combination of m off-resonance Raman transition and m control flelds, we can generate any \Dm) 
(or superpositions thereof). 


Preparation of many-body entangled states: error analysis. 


By using the effective hamiltonian ^eff = + He) ^ (with ^ being the hamiltonian projection into the DFS) under the 

conditions described in the previous Section, we see that the time of the operation to do a complete transfer of population from 
to by using an off-resonant transition, i.e., » ^tlsey is given by: 




(SM15) 


So far, we have considered the ideal situation; however, to estimate the fldelities in the preparation of these states, we need to 
analyze the errors that may occur within 4p ^ • The errors come from: 

• The spontaneously emitted photons from ) to other decay channels other than the waveguide, described by a Lindblad 
term: ^*(p) = Ln t {^geP^eg ~ P^ee) +H.C., that we embedded in a single decay rate F*. As we are using an off- 
resonant Raman transition, this source of error scales (for A ^ m > 1): 


= r 


mlClr 


4(yv„ + i)(A2 + (r*)2) 




m\Q,r\'^ 
4AA2 ’ 


(SM16) 


where we used that A^ ^ F*. 


• The other errors may appear due to photons emitted from the states out of the DFS, i.e., through the |F^^i)-like states. We 
can estimate the rate of these errors which are Anally given by (for A ^ m > 1): 


£ „ = (r+ArriD) 

A,S 

= (r*+A^riD) 


mlClr 


4A2(A2 + (r*+yvriD)2 


4(A2 + (r*+AriD)2) 


+ (r+yvriD) 


m|Qrp 

4;v(A2 + (r*+A?riD)2)) ’ 


(SM17) 
(SMI 8) 


Using these estimations, and with the following hierarchy for the parameters: A ^ m > 1 and Fid ^ A^ F* to simplify the 
expressions, we And the infldelity of the step m to be 

1 - F- = 4p^ (4"^ + f (^ + ^) ’ (SM19) 

which is optimized for A^ opt = VF*Fid yielding a scaling 1 / V^iD- Interestingly, the scaling does not depend on 

either the number of atoms. A, nor the number of excitations, m. 

So far, we have focused our discussion on the m-th step that goes from |D^-i) ^ |D^). As depicted in Fig. 2 of the main 
manuscript, by combining this process with a transition over the ancilla qubit that initializes it in \s)a, we can generate any 
arbitrary superposition of states (over the first A atoms) up to a mmax, namely I'Fd) = ^m\Dm), with a total of rnmax Raman 

transitions together with rnmax initialization gates on the ancilla. Neglecting the errors of the microwave transition over the 
ancilla, the total infldelity to generate I'Fd) is: 

1-Foptoc^. (SM20) 


Preparation of many-body entangled states: numerical analysis. 

In order to validate our scaling analysis, we study the preparation of two relevant sets of states without doing any approxima¬ 
tion and considering all the possible states (including super-radiant ones). The two sets of states are i) the general class of states 
\Dfn) and ii) the superpositions |d>^) = ^ (|Do) + |D^)). 
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First, it is interesting to realize that the symmetry conditions found from our analysis of the ideal situation tell us the relevant 
Hilbert space of the problem can be written in terms of the states: (8) | V^)a, where \ \1/)a = {|^)a, \s)a^ \^)a} is the Hilbert 

space of the ancilla atom, k can be restricted up to 2 photons (as higher excited states will be only weakly populated), and 
m = 0,1,..., mmax, where rnmax is the highest excitation that we want to achieve. Notice, that the Hilbert space does not depend 
directly on the total atom number A^, but only on rnmax- The number of atoms, N, only enters in the two-photon resonant condition 
that fixes . In this Hilbert space we can use a non-Hermitian evolution, where the effective Hamiltonian is determined by 
the action of Ffias + He plus the imaginary energies determined by the coupling to the waveguide through the Sge operator. The 
non-Hermitian Hamiltonian elements are given by: 


{^n,q\ ^ (0|A(F/las +^c)|Ti^,A:) ^ |V^)a — k[Ae — i{{N — m — k /2/2)]5k^q5fn,n^\f/^(j) + 


rklr 




-H.c.l 




^k,q^m,n^y/,s^^,e H.C.] 


• (SM21) 


A diagram of the relevant transitions in the complete Hilbert space is depicted in Fig. SMI (a). For the generation of individual 
Fock states, the pulse sequence can be easily deduced. One just needs to ensure a complete transfer of populations from 

|Tb,o) (8) |g)A ^ |Tb,o) (8) |*y)A ^ |T’i,o) (8) |g)A ^ |T’i,o) (8) |*y)A ^ 1^2,o) (8) |g)A-> |Tm,o) (8) |g)A, which can be done by fixing the 

time of interaction, t, to ^ = n) for the microwave (two-photon Raman) transitions. In Fig. 3(a) of the main text, 

we show the numerical fidelities obtaining when fixing the off-resonant transition to the optimal opt that we explored in the 
previous Section: the dots correspond to the numerical fidelities, whereas the solid lines depict the scaling cx 1 / ^/P^ and show 
how our general arguments give us the right scaling. For a more complicated state, such as the the pulse sequence can be 
calculated numerically. In Fig. 3(b), we show the optimal fidelities for generating these states up to 5 excitations, showing again 
how the 1 /\/Fid scaling of fidelities also holds for superpositions. 


Conditional preparation of many-body entangled states: using post-selection. 

In the error analysis we made in Section , we realized that some of the errors were coming from the small populations of 
superradiant states \Xg^'^) that emit quickly into the waveguide. Actually, in our atom-waveguide configuration one can think 
of using another atomic ensemble that acts as an efficient photonic absorber that maps the photonic excitation into an collective 
atomic one that afterwards can be detected via fiuorescence with a very high fidelity. Moreover, it is possible to use a more 
elaborate scheme as sketched in Ref. [25], where collective atomic excitations can be mapped to a single impurity atom, making 
fiuorescence detection much more efficient. As it is not the purpose of our manuscript to elaborate on conditional preparation, 
we just assume that we can perfectly detect all the photons emitted through the waveguide and study the scaling of the fidelities. 
Using these assumptions, we can cancel the errors in Eqs. SMI7 and SMI8 that are proportional to A^Fid, that in the limit 
A > m > 1 yield: 


~ ^ 4A^2(A2 + (r*+yvriD)2) ’ 

_ ^ m\Q.r\^ _ m\Q.r\^ _ 

xf’ ~ 4(a 2 + (r*+AriD)2) 4 a?(a 2 + (p +A?riD)2)) ‘ 


(SM22) 

(SM23) 


We can immediately see that the leading error comes from the first contribution of £ (m) - Taking the leading error only, we can 
then estimate the infidelity of the step m to be 




^(m) 


■ W ■ 


■ 


;rAr* / I Ae \ 

2 Ia;^ (AFiDpi’ 


(SM24) 


which is optimized for A^^opt = AFid yielding a scaling 1 — oc 1 /Fid. It is important to highlight that the optimal condition 
can not be realized as it implies A^ opt ^ Fid, which violates the conditions under which we derived our effective Hamiltonian. 
However, the linear improvement with 1 /Fid is still obtained even if we do not reach the optimal conditions. More details, on 
how to take advantage of the atom nanophotonic waveguide for conditional preparation will be presented elsewhere. [51] 
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ATOM-PHOTON MAPPING STARTING FROM AN INITIAL ATOMIC EXCITATIONS. 

Let us review first the general derivation for a Hamiltonian of the form H = Hsb , with 


+H.C. , (SM25) 

n,q 

where Hs (Hb) are the system (ID-bath) Hamiltonians, and ^sb is the interaction between them, where (aq) are the sys¬ 
tem (ID-bath) operators. Using the generalized input-output formalism [42], everything boils down to calculate the scattering 
amplitude 


A(r) = (0out|(5out|e-‘^'|Bin)|0in), (SM26) 

where |0in(out)) = denotes the system input (output) state at time t and |5in(out)) the input (output) state of the 

bath, which in our case is the electromagnetic field inside the waveguide. Let us particularize to our situation of interest, where 
we decay from an initial system state with m excitations, that is, ~ Tout{T) = 1 and ^in = 1 and the operators 

Moreover, we also linearize the waveguide dispersion relationship, i.e., (Oq ^ Vg{qsi)\q\, assume for simplicity that: 
gq ^ gq^ = g. With these considerations, the scattering amplitude for m excitations simplifies to [42, 52] 

A{,}(0 = x (vac|^0,, (^i)0,,(^2) ■ ■ ■ , (SM27) 

where {q} = {q \,..., is the set of relevant momenta of the m-photon state, where each of them run over the whole Brillouin 
Zone qi^B.Z. and ^ is the time ordering operator that guarantee that > ^'2 > • • • > A further simplification is obtained if 
we assume that we are within the Markov approximation, and use the fact that we work with an atomic configuration such that: 
q^Zn = 27tn. Then, 




n ^-iqaZn 


(SM28) 


With this approximation the time ordering simply ensures that the final (bosonic) state is symmetrized over all sets {q}. 
Notice, that the output photonic state associated to this scattering amplitude can be written as 

I^W) (SM29) 

where the sum over {q} = {qi,... ,qm} extends over all momenta. The state is normalized with the 1/m! factor as it 

cancels the m! terms that appears from the permutations of all the a^.’s and the m!-factor of the scattering amplitude normal¬ 
ization in the whole {^}-space, i.e., |A|^j(r)p = m! (notice the scattering amplitude is normalized to one only if: 

Y.qi>q2>--->qm \^{q}i^)\ ~ 

Therefore, it is enough to calculate the contribution of one time ordering, e.g., si > S 2 > ' > and then sum up to all the 

permutations of {q}. As was shown in the previous sections, the effective (non-hermitian) system Hamiltonian is 

^eff = (^di^ee — '^(XvD/2)SegSge. (SM30) 

Interestingly, our initial state \Sm) is an eigenstate of this effective Hamiltonian, and we can now calculate the action of the 
operator 


Sge{s)\S,n) = 


(SM31) 


and hence the correlator 


{vac\Sge{si)Sge{s2) 


= P[v^exp 

r=l 


-m 


'a 


riD(rA^r-(r-l)A^.-i)/2]^, 


(SM32) 


When doing the integral, one needs also to take care of the particular time ordering considered. In the case that we have chosen, 
i"! > • • • > 5-^, the integral can be rearranged Jq d^-^ dsm-i • • • dsi . The choice of the upper limit of integration to t is not 
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casual, as in each time integral we will obtain a term proportional to oc e disappear for times t ^ l/(A^riD) 

which are the ones that we are interested in. For example, the first integral: 


/- 

J S2 


d^-i exp 


- COa) -rijyNi/2]si 


exp 


-(0a)-riDM/2 
-(0a)-riDM/2)r ) 


exp 


(^exp (i(co^,-a^)-riDM/2)s2 


(SM33) 


-(0a)-riDM/2 


{-i{COq^ - (Oa) - rij)Ni/2)s2 


where the last approximation was done for r ^ 1 /(A^Fid). The second integral then reads: 

- fds2 - 

J Sl, 


ft exp + cOq^ - 2co^) - riD 2 A^ 2 / 2 ] ^ 2 ] 


i((0^l - (Oa) -riDA^l/2 
1 


(SM34) 


[i((0^i - (Oa) -rij^Ni/2][i{cOq^ + (0^2 “2(0a) -riD2A^2/2] 


exp 


i{C0q^ + (0^2 - 2^«a) - riD2A^2/2] ^3 


Iterating this integration and considering the permutations of the {q} due to the different time-orderings, we obtain the fol- 
lowing expression for the scattering amplitude 


\/rN~r e 




(SM35) 


for sufficiently large times t ^ I/A/^Fid, that is when the system state has completely decayed and all the excitations have 
been transferred to the bath state. Notice that the only dependence on t in this case enters through: e“‘Tr=i ^qrt^ which describes 
the center of mass motion of the wavepacket when going to the real space. In the low excitation regime, one can either do a 
Holstein-Primakoff approximation [45] or change Nm N, in the expression of Eq. SM35. In both cases we obtain: 




igVN 


- !(«,„ - «,)+r,„«75 = ■ 


(SM36) 


which represents a single mode wavepacket with spectral shape To emphasize the connection between the linear and 

non-linear scattering amplitudes we exemplify the results for the m = 2 photon wavepacket. The non-linear scattering amplitude 
is given in this case by 


^91,® (0 = 


+ 77 


+ ( 0^2 - 2ffla) -rir, 2 N 2/2 - O)^) - riDM/2 i{(Oq^ - ft)a) - riDM/2J 


= -g^^/2N2Nle~^^^r=l 


- r idM /2] - ft)a) - r idM /2] L 


i+riD7 


{N 2 -NI) 


{mq^+mq^-2ak)-TiD2N2/2\ ’ 

(SM37) 


where one can clearly see how if N 2 = Ni, then Aq^^q^{t) = Aq^q^{t). Moreover, d&Ni—N 2 = 1, then, it is direct to see that 
the correction due to the N 2 —N\ term is of the order 0{\/N). To generalize to higher photon numbers it is more convenient 
to directly calculate the overlap between the linear and non-linear approximations before doing the time integral. The reason is 
that one can formally integrate in [q] variables before the times Si, to obtain: 


E [>(%■ - f«a)(Sr - Sr)] = 


2L\g? 


5{Sr-Sr) =FiD5(^^-i'^). 


qi ^g(^a) 

With these 5’s the double time integral appearing when calculating the overlap is much simplified: 


(SM38) 






B I ^5,HP/ 


I 

{ 9 } 


= r 


ml 


ID 


rt rt m 

/ dSfn’” d^irTrV^^exp[-FiD(rA/^-(r-l)A/^_i)/2]]exp[-FiDrA^/2], 
Jo Js2 fLi 


(SM39) 
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where in the last equality we have used the fact that the contribution of m\ different time-orders will be the same. Then, the 
multi-time integral can be calculated iteratively, yielding to an overlap: 


1 - 




(m) , 

B,HP- 


= 1 - 2 " 


fJiN + Nr 


r=\ ^ ~ 


r/{2N) 


+ 0(mVA^^). 


(SM40) 


From the expression above, one can also check that is normalized by making Nr = N. For consistency, one can also check 

that each Cy^^n{^) in the linear expresion of is normalized independently: 


q 

-f 

Jo 


dq\g\^NL 


27t (a),-a)a)2 + (riDA^/2)2 


' d(0\g\^NL 


N 


0 Vg{(o^)n {(Og - (0^)^ + {TwN/2)^ 


L\g\^Nvg{( 0 ^) 2n 
It TidA^ 


= 1 


(SM41) 


IMPLEMENTATION DETAILS: PHOTONIC CRYSTAL WAVEGUIDES. 


A particularly promising system to implement our proposal is atom ID nanophotonic systems, in which first working proof- 
of-principles examples have been realized by using “alligator” photonic crystal waveguides [6, 10]. In these systems, the renor¬ 
malized spontaneous decay rate is given by: 


Fa “ 2A^ ’ 


(SM42) 


where Ug = c/vgh the group index, a = 3X^1 {In) the radiative cross-section. Am the effective mode area and ^ is an adimen- 
sional factor of cavity enhancement due the reflections at the end of the dielectric waveguide. Current SiN structures [6, 10], 
have Am ~ 0.2 /rm^, rig 10 and cavity enhancement ^ 5. There are several sources of errors in these systems: 


1. Spontaneous emission to other modes different from the chosen guided mode. Current structures show F* ^ F^, however, 
further design may result in further reduction of spontaneous emission, e.g., by making thicker dielectric structures F* ~ 
O.lFa [53]. Depending on the reduction of spontaneous emission F* = aFa, the Purcell factor with current designs can be 
PiD ^ 50/ Oi. 


2. Intrinsic losses of the material yield finite g-factors, which can be calculated as: 

Q=^ (SM43) 

2ni 

with n = nr — ini the refractive index of the material. The g-factor can be easily related to the attenuation of the intensity 
of the field traveling through the dielectric as follows: 




irop ' 


^ ^ g>^a 

AnniUg 2nng 


(SM44) 


with ^ 2 , = Xq!U r and where Lprop incorporates both material absorption via Hi and effect of reduced group velocity. We 
notice that g-factor also has contribution of scattering losses of the material due to imperfections, and therefore one 
must consider state-of-the-art values for doing estimations. For Cs atoms (Aq = 894 nm) and SiN structures {fiy = 2, 
g ^ 10^, = 10), this yields Lprop/Aa ^ 10"^- The main effect is that Jmn must be corrected by this attenuation length: 

Jmn ~ withz^„ = Zm—Zn- As Lprop/Aa IS Very large, in our situation of interest with Zn = nX^, 

the effect of the finite propagation can be treated as a perturbation to the collective Liouvillian given by: 


=Sfprop(p) = I^(l 


-e 1^"" 


l/Vop) (o-^.pa™ -pa^o-;,) +H.C. 


I 

n,m 


riD 


2 L, 


'prop 


((Tg>(T--p(T,Xe) +H.C.. (SM45) 


and introduces small corrections to the superradiant decay rate \Sm) and spontaneous emission rate of as long as 

the size of the atomic ensemble is NX^ ^ Lprop. 
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3. If one thinks of increasing Fid only through group velocity reduction another effect to take into account is retardation. 
The worst-case correction of this effect appears after doing a fast-resonant ;r/2 pulse to switch from |Z)^) to to do the 
atom-photon mapping. To observe superradiant behaviour in that case, it must be satisfied A^Fid < 2vg/(A^Aa) (assuming 
Aa/2 separation of atomic states), that with current state-of-art parameters [ 6 , 10], leads to N < 500. Notice, that in the 
preparation of superposition of |F>m), this critical number goes up to A^ < 10 ^, as the characteristic timescales do not show 
the collective enhancement. 

Furthermore, there are several ways of overcoming retardation in these set-ups: i) increase Fid not only by Vg but through 
cavity enhancement, e.g., by placing mirrors at the end of dielectric [ 10 ]; ii) more easily by reducing the characteristic 
timescale doing the atom-photon mapping from \Dm) to \Sm) off-resonantly by setting a finite 7 ^ 0. This reduces F^ 
(and F*) by a factor which relaxes retardation requirements, while keeping Fid constant. 

4. Moreover, a typical way of increasing group index is by using the regions of slow-light that appear close to ID band-gaps, 
where one can approximate the dispersion relationship by COq ^ COc+M^~ This dispersion will generate corrections 
with respect to the linear propagation of the wavepacket that must be kept small within the bandwidth mA^FiD- 
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